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A nonequilibrium, axisymmetric, Navier- Stokes flow solver with coupled radiation has been developed for 
use in the design of thermal protection systems for vehicles where radiation effects are important. The present 
method has been compared with an existing flow and radiation solver and with the Project Fire II experimental 
data. Good agreement has been obtained over the entire Fire II trajectory with the experimentally determined 
values of the stagnation radiation intensity in the 0.2 -6.2 eV range and with the total stagnation heating. The 
effects of a number of flow models are examined to determine which combination of physical models produces 
the best agreement with the experimental data. These models include radiation coupling, multitemperature 
thermal models, and finite rate chemistry. Finally, the computational efficiency of the present model is evaluated. 
The radiation properties model developed for this study is shown to offer significant computational savings 
compared to existing codes. 
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Nomenclature 

mass fraction for species s 

total energy 

vibrational energy 

enthalpy of species s 

surface intensity 

first-order reaction rate at wall 

Boltzmann constant 

mass of species s 

surface pressure 

source of radiation coupling 

source of translational -vibrational coupling 

heat flux 

total convective heat flux 

radiative heat flux 

translational heat flux component 

vibrational heat flux component 

characteristic temperature for dissociation 

wall temperature 

velocity component 

diffusion velocity component for species s 
chemical source term for species s 
direction component 
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y = fraction of species .v consumed at the surface 
p = mixture density 

Introduction 

A NUMBER of future NASA missions include scenarios 
where radiative heating may be important. For example, 
aerobraking scenarios use atmospheric drag in lieu of pro- 
pulsive braking for trajectory manipulation. For aerobraking, 
it is believed that the contribution of radiation to the surface 
heating may be 50% or more of the total. 1 Another class of 
missions is planetary probes. For planetary probes, the re- 
entry velocity is usually high enough to make radiation a 
significant or dominant heating mechanism. Therefore, to de- 
sign an appropriate thermal protection system (TPS), accurate 
predictions of the magnitude of the incident radiation and its 
interaction with the flowfield and the vehicle’s surface are 
needed. The major objective of this work is to develop a tool 
to aid TPS design for vehicles in which strong radiative effects 
are important. 

To develop a prediction tool for re-entry flows where ra- 
diation is important, four major areas are addressed. First, a 
flow solver is needed that incorporates the appropriate physics 
of the flow such as chemical and thermal nonequilibrium and 
viscosity. Second, radiation properties such as emission and 
absorption coefficients for the various radiation mechanisms 
must be determined. Third, a radiation transport algorithm 
is needed to describe the interaction between the emitted 
photons and the fluid flow. Finally, accurate boundary con- 
ditions are needed to describe the interaction of the fluid and 
the radiation with the vehicle’s surface. In this work, a two- 
dimensional or axisymmetric, Navier-Stokes flow solver that 
allows for chemical and thermal nonequilibrium is coupled 
with a detailed nonequilibrium radiation properties model 
that employs a tangent slab approximation for the radiation 
transport. 

A heating analysis for the design of a vehicle’s TPS gen- 
erally requires determining a total heat load for re-entry. 
Calculating a total heat load requires heating estimates at a 
number of trajectory points. This leads to an emphasis on 
computational efficiency. Towards this end, an effort has been 
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made to select the most efficient numerical algorithms. For 
the flow solver, a modified version of Candler’s nonequilib- 
rium reacting gas code 2 is used to generate Navier-Stokes 
flow solutions. In this study, this code will be referred to as 
GIANTS (Gauss-Siedel implicit aerothermodynamic Na- 
vier-Stokes code with thermochemical surface conditions). 
To determine the radiation properties, a new radiation pro- 
cess model has been developed from the LORAN code of 
Hartung.* The algorithms employed in LORAN have been 
restructured, rewritten, and modified to achieve maximum 
computational performance on a vector supercomputer. The 
resulting algorithm has been designated as the nonequilibrium 
optimized vectorizable radiation process model (NO VAR). 
The two codes, GIANTS and NOVAR, have been coupled 
together to calculate flow solutions about the Project Fire II 
vehicle 4 for a number of points along the vehicle’s trajectory. 
Numerical estimates of the heating and radiation intensity arc 
compared with the experimental data and previously gener- 
ated numerical solutions. Also, effects of thermal and chem- 
ical nonequilibrium are examined. Finally, the relative com- 
putational efficiency of the present algorithm is evaluated. 
Details of the methods are described in the following sections. 

Procedure 

Fluid Flow Governing Equations 

The appropriate fluid dynamic equations for a re-entry flow 
have been described in a number of sources. 2 s 7 As a start, 
a set of equations was chosen that includes most of the im- 
portant aspects of a re-entry flow in air. These equations 
include 1 1 species equations to determine the densities of N 2 , 
CL, NO, N, O, N t , O G NO + , N * , O + , and e , two total 
momentum equations, a total energy equation, and a separate 
energy equation for the vibrational energy. Solutions to this 
equation set for two-dimensional and axisymmetric flows are 
obtained using the Gauss-Seidel procedure of Candler and 
MacCormack. 2 H 

In this work, the thermal model w'as developed assuming 
that the translational temperature is equal to the rotational 
temperature and that one vibrational temperature is repre- 
sentative of all the vibrational modes of the gas mixture. The 
energy equations used to describe this model are as follows. 
The total energy equation is given by 

a(pg) al(pe + p)u\ d(u,r„ + q„ + q Vf ) 

at dx f te, 

[ "V/> ] 
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The vibrational energy equation is 
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where w s e„ is the vibrational energy lost or gained from chem- 
ical reactions. The coupling between the vibrational and trans- 
lational modes Q } r is modeled using a Landau— Teller for- 
mulation. The relaxation times for the vibrational modes are 
obtained from Millikan and White/' Radiation coupling enters 
the equations through the source term 0 rad in the total energy 
equation. In this study, the contribution of the electronic 
modes to the total energy and the radiative source of vibra- 
tional energy is neglected. For the radiation calculation, it is 


assumed that the electronic temperature is equal to the vi- 
brational temperature. This approximation greatly simplifies 
the computational effort. At the Fire II conditions where 
vibrational-electronic coupling and molecular chemistry ef- 
fects are important, this assumption may be a reasonable 
approximation. 

The chemical source terms are modeled using the Park 1993 
rates. 11 ’ However, for electron impact ionization, the rates 
deduced by Wilson 11 are employed. Also, the ionization rates 
are assumed to be a function of the translational temperature 
77 Ideally, these rates should be a function of the free electron 
temperature. Thus, this assumption should be applied with 
caution. In the reg ion of thermal nonequilibrium within the 
shock, T, t = \ / TT v is used as the controlling temperature for 
the dissociation rates. 

The mixture rules for the viscosity and translational and 
vibrational conductivity used in this work are based on the 
formulations of Gupta et alA 12 The cross section data by Yos 
in Refs. 5 and 12 has been reviewed and updated to improve 
the collision modeling. 14 The multicomponent diffusion coef- 
ficients arc treated using a bifurcation approximation. 14 |S 

Radiation Coupling and Transport 

In the Gauss-Seidel method, the chemical, inviscid and 
viscous terms are treated in an implicit manner. With the 
addition of radiation, a method of treatment, implicit or ex- 
plicit, must be selected for the radiation source term 0 rad . 
Gok^en treated the radiation source term in an implicit man- 
ner using the Gauss-Seidel method. 16 He found that an im- 
plicit treatment of the radiation source greatly increased the 
computational costs and reduced the stability of the overall 
method. For a detailed radiation model, the application of 
the method was limited to inviscid flows because of the in- 
creased computational cost. Thus, in the current work, the 
radiation is treated in a loosely coupled manner. 

The loosely coupled treatment of the radiation source has 
been described by Hartung et al. 4 The application of the loosely 
coupled method involves a number of steps. First, a converged 
flow solution is generated without radiation. The converged 
flow solution is then used to produce an initial estimate of 
the radiation field. Next, the flow r code is run for a certain 
number of iterations with the radiation field fixed. Then, the 
radiation field is updated using the new flow solution. This 
procedure continues until the flow and radiation solutions 
converge. 

The radiation transport is modeled using a one-dimensional 
tangent slab approximation. In the tangent slab model, ra- 
diation properties are assumed to vary only in one dimension. 
Perpendicular to the direction of variation, radiation prop- 
erties are assumed constant in planes that extend to infinity. 
For Fire II, tangent slab is a reasonable approximation in the 
stagnation region because of the highly blunted shape of the 
forebody. Most of the variation in flow properties occurs in 
the direction normal to the surface of the vehicle. 

Radiation Properties 

The radiation properties model employed in this study is a 
derivative of the Langley optimized radiative nonequilibrium 
code (LORAN). 17 1K LORAN generates absorption and emis- 
sion coefficients for the following radiative processes: atomic 
bound-bound, bound-free and free-free transitions for N 
and O and the N ; first negative, N. first and second positive, 
and Lyman-Birge-Hopfield, NO and y, and 0 : Schu- 
mann- Runge molecular band systems. The various atomic 
and molecular excited state populations are calculated using 
the quasi-steady-state (OSS) approximation of Park. 19 

When LORAN was developed, a major goal was to couple 
the method with a flowfield solver. The most sophisticated 
nonequilibrium radiation solver available is the NEQAIR code 
of Park. 21 ’ The NEQAIR code calculates the radiation prop- 
erties for the atomic and molecular processes listed in the 
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previous paragraph using a detailed line-by-line approach. 
This detailed approach was deemed too costly for flow cou- 
pling. Thus, a number of approximations were made to im- 
prove the computational efficiency of the method. 

The major differences between NEQAIR and LORAN are 
that the radiation properties for the molecular band systems 
are calculated using a “smeared band” approximation instead 
of the line-by-line approach used in NEQAIR, and an opti- 
mized spectral array is generated to reduce the total number 
of spectral points . 21 Thus, the accuracy of the method is con- 
sidered less than the detailed line-by-line approach of the 
NEQAIR model. However, the absorption and emission coef- 
ficients generated by the LORAN code have been shown to 
be comparable to those generated from the NEQAIR code 
of Park . 17 

Although the computational cost of the LORAN method 
is lower than that of NEQAIR and the feasibility of using the 
method in coupled calculations has been demonstrated, the 
computational cost to couple LORAN with a flow solver is 
still large. •' Initial timing tests on LORAN demonstrated that 
large portions of the code were nonvectorized. The low com- 
putational efficiency of the method was the major impetus 
for the development of the NOVAR algorithm. Timing com- 
parisons between NOVAR and LORAN are presented in the 
results section. 

The physical modeling of the various radiation mechanisms 
in LORAN and NOVAR are similar. However, to maximize 
the level of vectorization, the implementation of these models 
in the codes is quite different. A description of the algorithmic 
differences between LORAN and NOVAR are described in 
Ref. 15. 

One major difference between NOVAR and LORAN is 
the method used to calculate the populations of the excited 
molecular and atomic states. In LORAN, a QSS approxi- 
mation is employed to determine the populations of these 
states. In NOVAR, the QSS distribution is replaced by a 
Boltzmann distribution, which while theoretically not totally 
correct for the nonequilibrium region, may be adequate for 
the Fire II flight regime. 

In the course of developing the NOVAR algorithm, it was 
found that if an identical flow solution and excited state pop- 
ulation distribution were used as an input for both NOVAR 
and LORAN, then the differences between the calculated 
radiation properties were small. Figure 1 shows calculated 
absorption and emission coefficients from NOVAR and 
LORAN using equilibrium air at a pressure of 2 atm and a 
temperature of 10,000 K employing a Boltzmann distribution 
for the excited states. These conditions are representative of 


Equilibrium Air: P-2 atm, T- 10,000 K 



% Differer>ce for emission 



Fig. 1 Comparison of absorption and emission coefficients for equi« 
librium air conditions. 


the flow behind the shock for Fire II near peak heating. It 
can be seen in the figure that the percent difference between 
the calculated emission coefficients is on the order of 2 - 3 % 
or less. The radiation properties were generated at different 
spectral locations. Thus, some of the differences result from 
interpolating the emission coefficients onto an identical fre- 
quency distribution for the comparison. A similar result was 
obtained for the absorption coefficient. 

Boundary Conditions 

For this study, no-slip conditions are assumed. Thus, the 
velocity at the surface boundary is zero. The wall temperature 
is assumed to be constant, which is appropriate for the es- 
sentially “cold wall” of the Fire II experiment. Therefore, 

U s = U w = 0 T s = T w (3) 

where the subscripts s and w denote properties at the com- 
putational boundary and wall. 

For radiation, the wall is assumed to emit radiation as a 
blackbody at the specified wall temperature. The actual Fire 
II vehicle surface was nonblack. However, because of the low 
surface temperature, inaccuracies in this assumption are as- 
sumed small. In the freestream, the radiation flux is neglected, 
i.e., precursor effects are ignored. 

For chemistry, the wall is assumed to be fully catalytic to 
the ions and noncatalytic to the neutral species. Therefore, 
for the ions, the mass flux at the wall is related to the rate of 
diffusion towards the surface as follows 6 : 

(pD = ( Pv /C,)„. = h' vm , (4) 

where 17 denotes the direction normal to the surface. The 
reaction rate at the wall is given by 


~ y x V(kT w l2irm,) (5) 

where y s is the fraction of species s consumed at the surface, 
k is the Boltzmann constant, and m s is the mass of species s. 

For the Project Fire II experiment, the heat plug containing 
the thermocouples was constructed out of beryllium . 4 Beryl- 
lium is highly catalytic to reactions. Thus, a fully catalytic wall 
boundary condition has been assumed by many researchers 
when simulating the Fire II experiments . 2122 However, with 
the presence of atomic oxygen in the flow, it is likely that the 
surface of the beryllium was highly oxidized. In general, the 
catalytic potential of a material is related to its electrical prop- 
erties. Materials with a high electrical conductance such as 
metals usually act as catalysts for reactions. Conversely, ma- 
terials with a low electrical conductance, such as insulators, 
are poor catalysts . 22 Beryllium oxide has an extremely low 
electrical conductance and behaves as an insulator. Thus, in 
this work, the wall is assumed to be noncatalytic to reactions 
involving neutral species and y, for neutrals in Eq. (5) is set 
to zero. For the ions, the wall is assumed to be electrically 
uncharged. Thus, any ion reaching the surface is neutralized . 6 
Therefore, the wall is considered fully catalytic to ions and 

Yions 1 • 

Results 

The major objective of this study is to develop a coupled 
radiation solver for TPS design. To achieve this objective, a 
number of topics are discussed in this section. First, to assess 
the nonequilibrium flow modeling employed in GIANTS, an 
axisymmetric flow solution for the Project Fire II 1634-s con- 
ditions is generated and compared with results from the 
LAURA code of Gnoffo . 5 These flow calculations are used 
to generate uncoupled stagnation line radiation solutions us- 
ing NOVAR. Next, a baseline heating calculation is carried 
out over the Fire II trajectory using GIANTS/NOV AR. At 
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Table 1 Fire II trajectory conditions 


Time, 

s 

Altitude, 

km 

Velocity, 

km/s 

Density, 

kg/m A 


K 

r H „ 

K 

1634 

76.42 

11.36 

3.72 

X 

10 

S 

195 

615 

1636 

71.04 

11.31 

8.57 

X 

to 

s 

210 

810 

1637.5 

67.05 

1 1 .25 

1.47 

X 

10 

4 

228 

1030 

1639 

63.11 

11.14 

2.41 

X 

it) 

4 

242 

1325 

1640.5 

59.26 

10.97 

3.86 

X 

10 

4 

254 

1560 

1643 

53.04 

10.48 

7.80 

X 

10 

4 

276 

640 

1644 

50.67 

10.19 

1.02 

X 

10 

3 

281 

IKK) 

1645 

48.37 

9.83 

1.32 

X 

10 

3 

285 

1520 

1648 

42.14 

8.30 

3.00 

X 

10 


267 

1560 

1651 

37.19 

6.19 

6.05 

X 

10 

3 

253 

1060 


each trajectory point, uncoupled and coupled radiative fore- 
body flow solutions are generated. Then, numerically deter- 
mined values of the stagnation radiative intensity and the 
stagnation heat transfer are compared with the Fire 11 ex- 
perimental data and numerical results using LAURA/LOR AN. 
With a baseline established, the effects of a number of flow 
assumptions are examined. Radiative flow solutions are gen- 
erated with G I ANTS/N O V A R assuming thermal equilib- 
rium, and thermal and chemical equilibrium. Stagnation point 
results are compared with the experimental data and the base- 
line calculations. Finally, the computational efficiency of the 
present method is described. 

The flight conditions considered in this work are listed in 
Table l. 4 Progressing from 1634 to 1651 s, these conditions 
represent nonequilibrium to equilibrium flow conditions. Ex- 
perimentally measured values for the stagnation radiation in- 
tensity in the 0.2-6.2-eV range and total stagnation heating 
were extracted from plots contained in Rets. 4 and 24. In 
Ref. 25, the accuracy of the radiation data from the flight 
experiment was assessed and the data points were classified 
as prime or nonprime. Prime data points are considered to 
be the most accurate, i.e., where the flow was free of any 
contaminants. The Fire II flight experiment was composed of 
layers of beryllium and asbestos and the radiation sensor had 
a quartz window. For the nonprime points, it was believed 
the flow was contaminated or the quartz window w'as nearing 
its melting point. Thus, the radiation data at these points are 
not considered to be as accurate as at the prime data points. 

For all the GIANTS calculations in this work, grids with 
80 cells normal to the body and 37 cells along the body were 
employed. In past work, it was determined that this resolution 
was adequate . 2< ' The accuracy of this grid resolution including 
the normal spacing distribution is discussed in Ref. 26. 

Comparisons with LAURA 

For the first set of calculations, the nonequilibrium aspects 
of the present code are assessed. The amount of radiation 
emitted in the flowfield is directly proportional to the number 
of particles in each of the excited atomic and molecular states. 
The populations of these states are a function of the species 
densities and the vibrational and electronic temperature. In 
GIANTS, the electronic energy of the various species is ne- 
glected and the electronic temperature is assumed equal to 
the vibrational temperature for the radiation calculation. Thus, 
any of the flow mechanisms that affect the vibrational tem- 
perature can have a strong influence on the radiation calcu- 
lation. Two flow mechanisms that strongly influence the vi- 
brational temperature are the controlling temperatures used 
in dissociation and electron impact ionization. In GIANTS, 
the forward dissociation rates are assumed to be a function 
of the vibrational and translational temperatures and the elec- 
tron impact ionization reactions are assumed to be a function 
of the translational temperature. To illustrate the effects of 
these models on the flow and radiation calculations, a nu- 
merical solution for the 1634-s trajectory point is considered. 
For comparison purposes, a flow solution w r ithout coupled 


radiation from GIANTS is contrasted with a solution gen- 
erated from LAURA. In LAURA, the vibrational and elec- 
tronic modes are assumed to be a common energy pool char- 
acterized by one temperature T ir . The dissociation rates are 
assumed to be a function of the translational and vibrational- 
electronic temperatures and the electron impact ionization 
rates are assumed to be a function of the vibrational-elec- 
tronic temperature. 

In Fig. 2, temperature profiles along the stagnation stream- 
line are compared from GIANTS and LAURA. It can be 
seen in the figure that LAURA predicts a much higher trans- 
lational temperature in the nonequilibrium portion of the flow. 
Both methods predict a similar equilibrium temperature be- 
hind the shock. The major reason for the large differences 
between the nonequilibrium temperatures is the treatment of 
electron impact ionization. In LAURA, using T vr as the con- 
trolling temperature in the electron impact ionization rates 
lowers the amount of ionization, which leaves more energy 
in the translational modes. The vibrational-electron-elec- 
tronic temperature is increased via translational -vibrational 
and translational heavy particle-free electron coupling, which 
leads to an overshoot of the vibrational-electron-electronic 
temperature behind the shock. In (HANTS, using T as the 
controlling temperature in the electron impact ionization rates 
increases the amount of ionization, which low r ers the trans- 
lational temperature. The effect of using different rate con- 
trolling temperatures on the species densities is illustrated in 
Fig. 3, which shows the nitrogen species, N : , N, N>, and N 
along the stagnation streamline. The faster rate of N f ion 



Fig. 2 Temperatures along the stagnation line for 1634-s, axisym- 
metric flow. 



Fig. 3 Nitrogen species along the stagnation line for 1634 s. 
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production is clearly shown for GIANTS. Also, it can be seen 
that LAURA predicts a much higher N ^ concentration in 
the nonequilibrium portion of the flow. The larger N % density 
for LAURA results from the high nonequilibrium transla- 
tional temperature, which promotes charge exchange reac- 
tions that produce N 2 . In the region of thermal equilibrium 
behind the shock, both methods predict similar nitrogen spe- 
cies concentrations. Near the surface differences in the species 
densities can be attributed to the wall catalysis boundary con- 
ditions employed in the methods. A noncatalytic wall for the 
neutrals is assumed in GIANTS, while a fully catalytic wall 
for all species is assumed in LAURA. 

The effects of the nonequilibrium flow modeling on the 
radiation calculations are illustrated in Figs. 4 and 5. Uncou- 
pled stagnation line radiation solutions are generated from 
the GIANTS and LAURA flow solutions using NOVAR. 
Figures 4 and 5 show the radiative flux q r away from the wall 
and toward the wall along the stagnation line. In the tangent 
slab approach, the net radiative flux is the difference between 
the flux towards and the flux away from the wall. First, it can 
be seen that the radiative flux generally follows T v and T ve . 
In Fig. 4, the radiative flux away from the body rises as the 
vibrational temperature increases and levels off at the shock 
as the temperature drops. The overshoot in T ve for the LAURA 
calculation leads to an overshoot in the radiative flux. In Fig. 
5, the radiative flux toward the body quickly increases as T v 
and T ve rise. Nearest the wall, the radiative flux drops slightly 
with temperature as the result of absorption in the thermal 
boundary layer. 



Fig. 4 Radiative flux away from the wall along the stagnation line 
for 1634 s. 



Fig. 5 Radiative flux toward the wall along the stagnation line for 
1634 s. 


Table 2 Comparison of radiative intensities for Fire II 1634 s 



Avail ' 

0.2-16.5 eV, 
W/cirF-sr 

/wall 1 

0.2-6. 2 eV, 
W/cm 2 -sr 

^ wall » 

W/cm 2 

Experiment 

— 

1.3 ± 20% 



LAURA/NOVAR 

12.8 

4.2 

67.9 

GIANTS/NOVAR 

4.7 

1.4 

24.8 



Fig. 6 Comparison of total stagnation heating down the Fire II tra- 
jectory. 


The higher radiative flux for LAURA results mostly from 
two radiation mechanisms, atomic line and molecular NJ. 
The higher nonequilibrium T ve profile increases the popula- 
tions of the higher electronic states of N and O, which results 
in more atomic line radiation. A secondary effect is molecular 
radiation from NJ, which is a strong radiator. As shown in 
Fig. 3, the density of N 2 f is much larger for LAURA than 
for GIANTS in the shock region. Thus, the higher Nj density 
and larger T ve profile increase the molecular radiative flux. 

Numerical and experimental surface radiation intensities 
and heating values for the two uncoupled flow solutions are 
summarized in Table 2. In the experiment, the stagnation 
radiation intensity in a frequency range between 0.2-6. 2 eV 
was measured. The experimental value for 1634 s was 1.3 W / 
cm 2 with an error range of ±20%. The uncoupled surface 
intensities generated from the LAURA solution overpredict 
the experimental surface intensity in the sensor range by a 
factor of about three. The intensity value generated from 
GIANTS is in good agreement with experimental data. The 
LAURA calculation also produces a much higher fraction of 
uv radiation. This radiation comes primarily from atomic lines 
in the shock. Although some of this radiation is self-absorbed, 
a significant portion is transported in the wings of the atomic 
lines, 27 which results in a higher radiative surface heating 
value. 

Trajectory Calculations 

With a baseline for comparison established between the 
GIANTS and LAURA flow codes, coupled radiative calcu- 
lations using GIANTS/NOVAR are carried out over the en- 
tire Fire II trajectory. Numerically determined values for the 
uncoupled and coupled total stagnation heating values and 
radiative intensities are compared with the Fire II experi- 
mental data and with numerical values obtained using LAURA/ 
LORAN. The LAURA/LORAN data were extracted from 
plots published in Ref. 21. Figure 6 is a comparison of the 
experimental total stagnation heating (convective and radia- 
tive) along the Fire II trajectory with uncoupled and coupled 
radiation solutions using GIANTS/NOVAR and LAURA/ 
LORAN. The total stagnation heating is composed of the 
convective contribution and the radiation absorbed into the 


OI.YNK K FT A1 


591 


surface accounting for the absorption properties of beryllium. 

It can be seen that there are large differences between the 
predicted heating values from the methods. In the early por- 
tions of the trajectory both methods underpredict the total 
heating. It appears from 1634-1637.5 s that both methods are 
in good agreement. However, a fully catalytic wall boundary 
condition is assumed for LAURA, while a noncatalytic 
boundary condition for the neutrals is assumed for GIANTS. 
If identical boundary conditions were used for both methods, 
then the heating would be much higher for the GIANTS 
solution in the early part of the trajectory and the discrep- 
ancies between the methods would be larger. If the surface 
of the vehicle was only partially oxidized, then assuming a 
small finite catalytic boundary condition for the neutrals in 
GIANTS would improve the agreement with the early part 
of the flight experiment. Near peak heating, the coupled 
GIANTS solution tracks the experimental data fairly well. 
The LAURA results underpredict the data. It can be seen 
for the GIANTS results that coupling the radiation with the 
flow improves the agreement with the experimental data by 
lowering the amount of radiative heating. This effect is illus- 
trated in Fig. 7. 

Figure 7 shows a comparison of the measured and predicted 
stagnation intensity in the 0.2-6.2-eV range. The solid circles 
denote prime data points where the experimental data were 
considered to be most accurate. It can be seen in the figure 
that predicted intensity values from GIANTS/NOVAR are 
in good agreement with the experimental data near the prime 
data periods. In the 1 643— 1648-s time period, the GIANTS 
solution overpredicts the radiation intensity. However, the 
accuracy of the radiation data in this period is suspect. Fur- 
ther, the agreement of GI ANTS with the total heating in this 
period was good. The LAURA/LORAN combination over- 
predicts the radiation at 1634 s and underpredicts the radiation 
near its peak at 1643 s. However, agreement with the exper- 
imental data is good in the 1643- 1648-s time period over- 
predicted by the present method. It can be seen for both 
methods that radiation coupling lowers the amount of radia- 
tion. As expected, the coupling effects are largest where the 
radiation intensity is largest. Discrepancies between the 
GIANTS/NOVAR and LAURA/LORAN solutions could be 
the result of the different nonequilibrium models employed 
in both LAURA and GIANTS, as previously demonstrated. 
Another factor could be the OSS distribution that is employed 
in the LORAN code. Some effects of the OSS approximation 
on the radiation calculation are examined in Ref. 15. 

Figure 8 shows a comparison of the convective heating pre- 
dicted by GIANTS and LAURA along the flight trajectory. 
Experimental data w r ere not available for the convective heat- 
ing. It can be seen that large differences exist in the convective 



Fig. 7 Comparison of stagnation intensities down the Fire II trajec- 
tory. 



Fig. 8 Comparison of convective heating down the Fire II trajec- 
tory. 


heating predicted by both methods. Some of the discrepancy 
can be attributed to differences in the diffusion models em- 
ployed in GIANTS and LAURA. However, these differences 
are magnified by the fact that different catalytic boundary 
conditions were used for both LAURA and GIANTS and 
that radiation coupling had only a small effect on the con- 
vective heating. The largest differences occur at peak heating. 


Nonequitibrium Modeling Effects 

In the previous set of calculations a stagnation heating anal- 
ysis was generated for the Fire II trajectory using a set of 
nonequilibrium models that included a two-temperature ther- 
mal model and finite rate chemistry. Using these models in 
GIANTS/NOVAR, a favorable comparison was obtained with 
the Fire II experimental data. In the next set of calculations, 
simpler models are used to describe the nonequilibrium as- 
pects of the flow. Coupled radiative calculations are generated 
using GIANTS/NOVAR assuming a one-temperature ther- 
mal model with finite-rate chemistry and assuming a one- 
temperature thermal model and chemical equilibrium. These 
flow solutions were carried out to determine if less complex 
flow modeling could be used over all or part of the Fire II 
trajectory and still yield satisfactory agreement with the ex- 
perimental data. 

In GIANTS, the effects of a one-temperature thermal model 
are simulated by decreasing the relaxation time between the 
vibrational and translational modes by a factor of 3000. To 
simulate thermal and chemical equilibrium, the vibrational 
relaxation times are decreased by a factor of 3000 and the 
chemical rates are increased by a factor of 1000. Flow solu- 
tions using a simulated one-temperature thermal model and 
simulated equilibrium are generated with and without radia- 
tion for three trajectory points, 1634, 1643, and 1651 s. These 
trajectory points represent nonequilibrium, near-equilibrium, 
and equilibrium flow conditions. 

Figure 9 is a comparison of total stagnation heating for the 
one-temperature, equilibrium, and previous nonequilibrium 
calculations with the Fire II experimental data. The total heat- 
ing predicted by the one-temperature and equilibrium formu- 
lations is higher than the nonequilibrium formulation and the 
experimental data at each of the three calculated trajectory 
points. A one-temperature thermal model increases the total 
heating because the translational rotational — vibrational 
temperature overshoots the postshock temperature. The in- 
creased temperature leads to an increase in the radiative heat- 
ing. The equilibrium model increases the heating because the 
faster chemical rates in the cold boundary layer increase the 
rate of atomic recombination, which raises the translational 
temperature. The increased translational temperature leads 
to higher convective heating. These effects are illustrated in 
Figs. 10-13. 
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Fig. 9 Effects of thermal and chemical equilibrium on the stagnation 
heating along the Fire II trajectory. 
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Fig. 10 Vibrational and total temperatures in the shock region for 
the 1643-s conditions. 
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Fig. 1 1 Effects of thermal and chemical equilibrium on the absorbed 
radiative heating. 


Figure 10 is a plot of the vibrational temperature and the 
translational-rotational-vibrational temperature in the shock 
region for the 1643-s conditions. It can be seen that a one- 
temperature model overshoots the postshock temperature that 
all three formulations approach. As previously described in 
Figs. 3-5, an overshoot in the vibrational temperature can 
produce a large increase in the radiative flux. The overshoot 
for the one-temperature model is the result of chemical non- 


equilibrium in the shock region. This result was also observed 
for the 1634-s conditions. At the 165 1-s conditions, the flow 
is near chemical equilibrium in the shock and all three models 
yield similar temperature profiles. In contrast, the equilibrium 
temperature is much closer to the nonequilibrium vibrational 
temperature that was observed at all three trajectory condi- 
tions. 

Figure 11 shows the radiative heating absorbed by the be- 
ryllium surface for all three formulations. As expected, a one- 
temperature thermal model results in a significantly larger 
absorbed radiative flux at the 1634- and 1643-s conditions. 
The uncoupled radiative flux was so large at the 1634 s con- 
ditions that a coupled radiation solution was not generated. 
The absorbed radiation for the equilibrium formulation is 
similar to the nonequilibrium formulation at each of the cal- 
culated points. At the 1 65 1-s conditions, the flow is in chem- 
ical and thermal equilibrium and all of the formulations pro- 
duce similar results. 

Figure 12 shows a comparison of convective heating along 
the trajectory for the three formulations. The convective heat- 
ing for the one-temperature thermal model is similar to the 
nonequilibrium formulation at each of the calculated points, 
while the equilibrium formulation is higher. The equilibrium 
formulation leads to higher convective heating because the 
accelerated reaction rates result in more atomic recombina- 
tion in the boundary layer. The release of energy in the bound- 
ary layer raises the translational temperature and the con- 
vective heating. The increased recombination in the boundary 
layer is illustrated in Fig. 13, which shows a log- log plot of 
the mass fractions of N 2 , 0 2 , N, and O near the stagnation 



Fig. 12 Effects of thermal and chemical equilibrium on the convective 
heating along the Fire II trajectory. 
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Table 3 Comparison of computational 
efficiency on a Cray*C90 
(50 radiation points) 



Total CPU 



time, s 

Mflops 

LORAN 

53.43 

61.6 

NOVAR 

8.37 

332.8 



Fig. 14 Convergence history for the Fire II-1639-s calculation. 


point at the 1643-s conditions. It can be seen that the flow is 
close to equilibrium about 1 cm from the wall. Then, the 
nonequilibrium and equilibrium formulations diverge as the 
cold wall is approached. Even at the 1 65 1 -s conditions, this 
effect is observed because there is still atomic oxygen present 
in the boundary layer. More details about these calculations 
can be found in Ref. 15. 

Computational Efficiency 

One of the objectives of this study was to develop a tool 
for the prediction of radiating flowfields that is both accurate 
and numerically efficient. At the start of this study, timing 
estimates of the LORAN code indicated that its performance 
was poor on a vector machine. Since calculating the radiative 
properties is a large fraction of the total computational cost, 
an effort was made to improve the computational efficiency 
of the radiation algorithm, which resulted in the NO VAR 
code. Table 3 shows the CPU time required to calculate the 
radiation properties and the radiation transport for 50 points 
along the stagnation line at the 1634-s conditions. Timing 
comparisons were performed on a Cray-C90. The present 
method requires approximately 8 s to perform this calculation. 
To perform a similar calculation with LORAN required about 
53 s. Thus, over a factor of 6 improvement in the computa- 
tional efficiency has been obtained with the radiation solver 
developed for this work, NOVAR. The reduced computa- 
tional time for NOVAR is the result of extensive vectorization 
and optimization for the Cray-C90 and algorithm modifica- 
tions to make the method more efficient. 

Figure 14 shows a convergence history of the 1639-s cal- 
culation for the GIANTS/NO VAR combination. For the first 
1200 iterations, GIANTS is run without radiation, then the 
first radiation solution is generated using NOVAR. GIANTS 
is then run for 200 flow iterations with the radiation solution 
fixed. The calculation of a new radiation field and 200 flow 
iterations constitutes one radiation iteration. Each of the spikes 
in the residual represent the calculation of a new radiation 
field. It can be seen that the radiation field is close to con- 
verged by the fifth radiation iteration because the calculation 
of a new radiation solution has only a small effect on the 
overall residual. The CFL number after the 450th iteration 
was 400. The CFL number for the flow solver was not lowered 


when radiation was included in the calculation In general, it 
was found that the stability of the flow code was unaffected 
by radiation coupling for the cases considered in this study. 
The total CPU time used for this calculation (uncoupled and 
coupled) was about 2100 s on a Cray-C90. The calculation of 
the radiation properties required about 20 °/r of the total CPU 
time. 

Concluding Remarks 

A nonequilibrium, axisymmetric, Navier-Stokes flow 
solver with coupled radiation, GIANTS/NO VAR, has been 
developed for use in the design of thermal protection systems 
for vehicles w here radiation effects are important. The present 
method was compared with an existing flow and radiation 
solver and with the Project Fire II experimental data. Ex- 
cellent agreement was obtained over the entire Hre II tra- 
jectory with the experimentally determined values of the stag- 
nation radiation intensity in the 0.2-6. 2 eV range and the 
total stagnation heating. 

In comparing with the experimental data, the sensitivity of 
the numerical simulation to various flow models was exam- 
ined. From this sensitivity study a number of conclusions can 
be made about what flow models are needed to accurately 
simulate the Fire II experiment. A minimum combination of 
a two-temperature thermal model consisting of a separate 
translation -rotational temperature and a vibrational temper- 
ature was needed to obtain a good agreement with the ex- 
perimental radiation intensity. It was found that the popu- 
lations of the various vibrational and electronic states could 
be characterized by a Boltzmann distribution at the vibra- 
tional temperature. Thus, including the electronic energy and 
solving for a separate electronic temperature was unneces- 
sary. However, employing a one-temperature thermal model 
was found to be inadequate for simulating the experiment. 
For trajectory conditions with large thermal nonequilibrium, 
the one-temperature model overpredicted the experimental 
stagnation radiation intensity. 

To simulate chemical effects, finite rate chemistry with elec- 
tron impact ionization characterized by the translational tem- 
perature provided the best agreement with the experimental 
data over the entire trajectory. Whether or not this result is 
fortuitous or due to correct physical modeling is not known 
at this time. A flow calculation simulating equilibrium con- 
ditions was found to overpredict the experimental stagnation 
heating. The increased heating resulted from excessive atomic 
recombination in the boundary layer that increased the con- 
vective heating. 

The computational efficiency of the radiation properties 
model developed in this work, NOVAR. has been signifi- 
cantly improved over LORAN. The present model was found 
to be about six times faster than LORAN on a Cray-C90. 
Finally, the stability of the Navier-Stokes solver, GIANTS 
was not affected when loosely coupled with radiation. 
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